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Abstract 

The dinitrides CN2, SiN2, and GeN2 in assumed pyrite-type structures are studied 
by means of density functional theory using both ultrasoft pseudopotentials and 
the augmented spherical wave (ASW) method. The former two materials constitute 
the large-x limit of the broader class of CNx and SiNx compounds, which are well 
known for their interesting mechanical and electronic properties. For CN2 a large 
bulk modulus Bq of 405 GPa was determined . While SiN2 is found to be a wide 
band gap compound, the calculated gaps of CN2 and GeN2 are considerably smaller. 
The trends in structural and electronic properties as e.g. bond lengths, band gaps 
and covalency are well understood in terms of the interplay of different types of 
bonding. 
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1 Introduction 



In recent years nitrogen lias turned out to play a key role in todays most ex- 
citing technological applications of materials science. Important examples are 
the developement of III-V nitrides (AIN, GaN, InN) for light emitting diodes 
(LED) and lasers [1,2], the investigation of Si3N4 and SiOxNy as promising 
candidates for ultrathin high dielectric materials substituting Si02 in the race 
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for ever shrinking field effect transistor (FET) gate lengths [3], and, finally, 
new ultra hard materials with the proposal of a model carbon nitride C3N4 as 
the hardest known material [4]. However, all attempts to grow thin layers of 
nitrogen rich stoichiometrics failed because of the easy release of nitrogen due 
to the high stability of molecular N2 [5]. This major drawback led to propose 
nitrogen poor stoichiometrics such as C11N4 [6]. Nevertheless, the process of 
N2 release from CNx has not yet been investigated in detail and only recently a 
preliminary approach based on a model structure CN2 was presented [7]. How- 
ever, our knowledge about this class of nitride systems is far from exhaustive. 
To some extent this is due to the fact that the search for compounds with 
new materials properties is not straightforward and can be experimentally 
very demanding. In this situation, state-of-the-art computer simulations built 
on the principles of quantum mechanics as embodied in density functional 
theory (DFT) are very helpful. In particular, they even allow to investigate 
compounds, which are not yet synthesized. 

Aiming at a deeper understanding of ultra hard nitrides we report in the 
present work on investigations of the lightest IV-V compounds CN2, SiN2 and 
GeN2 as crystallizing in an assumed cubic pyrite structure [7]. This choice 
of hypothetical structure is guided by the observation that it allows to build 
up C-N bonds in a three-dimensional network while still satisfying the strong 
N-N bonding. At the same time, the simplicity of the pyrite structure offers 
the advantage that the different types of bondings can be easily addressed. 
Finally, this structure is also adopted by other isoelectronic AX2 compounds, 
as e.g. SiP2 [8,9]. Further experimental evidence comes from N- or C-doped 
SiP2 and related compounds. In the pyrite structure, the N-N entities occupy 
octahedral holes formed by the C, Si and Ge fee sublattices. The latter are 
well known from the cubic structures of diamond, SiC, SiGe, and TiC. As 
estimates for amorphous SiNx [10] and SiN2 clusters [11] show, CN2 and SiN2 
represent the large-x limit of those ultra hard and wide electronic band gap 
CNx and SiNx materials, which contain N-N-bonds. 

From a chemical point of view, the transition from SiP2 to SiN2, CN2, and 
GeN2 involves several fundamental issues. These are the behaviour of the light 
group IV compounds C, Si and Ge in octahedral environments as well as that 
of nitrogen in the center of tetrahedra formed by one N and three C (Si/Ge) 
neighbours. The resulting description on the related bonding schemes should 
be found in between Paulings covalent bonding [12] and ionic concepts like 
A^+[N2]^~ with the anions isoelectronic to [S2]^~ in pyrite FeS2. While [N2]^~ 
was proposed to appear in diacenides [13], covalent N-N-bonds have been pro- 
posed for clustered III-V semiconductor nitrides [14]. 

In this work we concentrate, using DFT-based methods, on finding the opti- 
mal crystal structure parameters of pyrite CN2, SiN2 and GeN2. In addition, 
we aim at relating the crystal structure stability to the electronic properties 
as well as to the chemical bonding. 
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2 Computational methods 



The electronic structure calculations were performed within the framework 
of density functional theory and the local density approximation using two 
complementary approaches. First, full geometry relaxations including opti- 
mization of the cell volume and the atomic positions were performed with the 
help of an ultrasoft pseudopotential scheme as implemented in the VASP code 
[15,16,17,18]. Calculations were converged at an energy cutoff of 434.82 eV for 
the plane- wave basis set. The Brillouin zone integration was performed using 
11 to 76 irreducible k points and the tetrahedron method [19]. 
In a second step, the electronic properties of all three compounds were de- 
termined using the all-electron augmented spherical wave (ASW) method 
[20,21]. In order to represent the correct shape of the crystal potential in 
the large voids, additional augmentation spheres were inserted into the open 
pyrite structure. Optimal augmentation sphere positions as well as radii of 
all spheres were automatically generated by the sphere geometry optimization 
(SCO) algorithm [22]. The basis sets comprised s, and 3d states of the group 
IV element, N 2s, 2p states as well as states of the additional augmentation 
spheres. The Brillouin zone integration was done using an increasing number 
of k points within the irreducible wedge ranging from 11 to 1135 points, again 
to ensure convergence of the results with respect to the k-space grid. 



3 Results and Discussion 

3.1 Crystal structure optimization 

The pyrite structure AX2 is best described in terms of the rocksalt structure 
with the A atoms and X2 occupying the two sublattices such that cubic sym- 
metry with space group Pa3 is preserved [23]. As a result, the A atoms are 
found at the center of octahedra, while three A atoms and one X atom form 
distorted tetrahedra about the X sites, see Fig. 1. The pyrite structure is fully 
specified by two parameters, namely, the lattice constant and the positional 
parameter x of the X atoms. 

To start with, we performed a structural optimization for SiP2, where the 
theoretical results could be checked against experimental data [8], see Tab. 1. 
As expected, the lattice constant is slightly underestimated by the LDA cal- 
culations. In contrast, the internal parameter xp is almost exactly reproduced 
and so is the fraction 4^^, with the P-P distance being shorter than the Si-P 
bond. 

On going to the nitride we observe a drastic volume decrease with a calcu- 
lated lattice constant of ~ 4.42 A close to the value of 4.348 A measured for 
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Fig. 1. Pyrite structure. Group IV element atoms are found at the centers of the 
cube and the edges, respectively. 
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Table 1 

Crystallographic data of the optimized structures. Lengths in A and energies in 
eV/atom. Experimental data taken from Ref. [8]. 

SiC [24]. In addition, the distance fraction has considerably decreased. 

While the N-N bond length amounts to 1.454 A, the Si-N distance of 1.886 A 
is very similar to the corresponding value of 1.889 A observed for Si3N4 and 
is indicative of a single bond. For GeN2, our calculations yield a bigger cell, 
albeit, with a still shorter N-N distance of 1.428 A. 

Finally, optimizing the pyrite structure for CN2, we obtain an even smaller 
lattice constant of 3.982 A. Together with the calculated internal parameter 
of xn = 0.403 this leads to an N-N distance of 1.34 A. While still being much 
larger than in molecular N2 (1.09 A), this value is considerably smaller than 
in SiN2 and GeN2 strengthening the double bond. At the same time the C-N 
distance of 1.695 A is significantly larger than the 1.47 A measured for the 
C-N single bonds of C3N4 indicating an even lower order. Applying the usual 
estimates for the hardness of CNx compounds [25,26] we calculate a bulk mod- 
ulus of Bq = 405 GPa from a Birch fit [27]. Since this is harder than cubic 
BN, CN2 in the pyrite type structure may be regarded as a high-density ul- 
tra hard material. Since the energy values given in Tab. 1 are relative to the 
atomic energies and thus not indicative for the stability of the compounds, we 
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performed additional supercell calculations for atomic C (Si, Ge) as well as 
molecular N2 and obtained a heat of formation of -4 eV/atom. 

3.2 Electronic structure 

The band structure of CN2 along selected high symmetry lines within the first 
Brillouin zone of the simple cubic lattice [23] as well as the partial densities of 
states (DOS) of all three compounds are shown in Figs. 2 and 3. In these figures 




Fig. 2. Electronic structure and partial densities of states (DOS) of CN2. 
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Fig. 3. Partial densities of states (DOS) of SiN2 and GeN2. 

four groups of bands are identified. The two low-energy groups of bands each 
comprise four single bands. Bands are most easily counted along the line M-R, 
where they are fourfold degenerate. According to the partial densities of states 
these bands trace back mainly to the N 2s states, which form a bonding and 
antibonding subset of several eV width. Whereas the large splitting between 
these two groups is due to the short N-N distance the intrinsic band widths of 
both groups mainly result from the dispersion across the face-centered cubic 
(fee) lattice formed by the nitrogen pairs. This is concluded from the level 
sequence at the F point into a lower nondegenerate and an upper threefold 
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degenerate state. The latter results from folding the fee bands into the simple 
cubic Brillouin zone (for the antibonding group of bands the sequence of levels 
is reversed). Further support comes from the striking similarity of the density 
of states especially of the bonding bands to that of a band calculated within 
the tight-binding approximation for an fee lattice. Band widths are largest for 
CN2, which has the smallest lattice constant. 

The electronic structure of the valence and conduction band are influenced by 
different types of bondings. Both groups comprise 20 bands, which derive from 
hybridization of all orbitals except for the N 2s states. Note that we have used 
a rotated coordinate system to separate the N 2p partial densities of states 
into contributions from orbitals mediating a- and vr-type bonding within the 
N2 pairs. The calculated optical band gaps are 0.9, 5.5, and 1.4 eV for CN2, 
SiN2 and GeN2, respectively. 

Common to all three nitrides is the strong bonding between group IV element 
p states and those of nitrogen, which leads to bonding and antibonding bands 
in the energy regions [-12:-6] and [5:9] eV for CN2, [-10:-4] and [10:15] eV for 
SiN2, and [-9:0] and [5:18] eV for GeN2, respectively. Since the nitrogen atoms 
are coordinated by three group IV atoms this bonding is rather isotropic and, 
hence, involves both the 2p^ and 2^,^ states in the same manner with the ratio 
of 1:2 of their partial DOS reflecting their respective weight. Except for the 
carbon compound the unoccupied antibonding bands are dominated by the 
group IV element p states, whereas the N 2p orbitals show a larger contribu- 
tion to the bonding bands. 

Differences between the three compounds, which cause, in particular, the 
spread in calculated optical band gaps, can be traced back mainly to the 
different influence of two additional types of bonding. These are, respectively, 
the bonding between the group IV element s and the N 2sp states as well as 
the strong a- and vr-type overlap within the nitrogen pairs. The former leads 
to considerable contributions of the C/Si/Ge s states to the second group 
of bands as well as to the lower edges of the valence and conduction bands. 
Across the series, s-p hybridization of the group IV element decreases and, on 
going from C to Ge, the bonding and antibonding bands involving the s states 
are lowered in energy relative to the bands derived from the group IV element 
p states. This effect is already visible for SiN2 and becomes most obvious for 
GeN2, where it leads to the substantial decrease of the optical band gap. 
To the contrary, the a- and yr-type overlap within the N2 pairs causes splitting 
of the 2pu and 2]?^ states into bonding and antibonding bands. In the partial 
DOS, those N 2p states, which are subject to overlap within the pairs, can be 
easily identified by the smaller C/Si/Ge contributions in the respective energy 
range. However, due to the large N-N distance in both SiN2 and GeN2, this 
bonding/antibonding splitting is rather weak in these two nitrides. This is con- 
trasted by the situation in CN2, where the close coupling of the nitrogen atoms 
leads to strong splitting into bonding and antibonding bands showing up in 
the energy intervals from [-6:0] and [8:12] eV, respectively. To conclude, the N 
2p states are subject to two different types of bonding, namely, that inside the 
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pairs and that with the group IV element sp states. The bonding and anti- 
bonding states resuhing from both types of chemical bonds end up in different 
energy regions. In particular, for CN2 the bonding bands resulting from N-N 
2p overlap are higher in energy than the bands growing out of C/Si/Ge sp-N 
2p overlap. As a consequence, the optical band gap is considerably reduced as 
compared to SiN2. 



4 Conclusions 



First principles calculations for the dinitrides CN2, SiN2 and GeN2 in assumed 
pyrite-type structures result in stable crystals. A bulk modulus of Bq = 405 eV 
is calculated for CN2, which thus can be regarded as an ultra hard material. 
The electronic structures of all three nitrides are governed by three different 
types of bonding. Overlap of the C/Si/Ge s and p orbitals with the N 2p states 
lead to the formation of bonding and antibonding bands, which give rise to 
the valence and conduction groups of bands, respectively. In both groups, the 
C/Si/Ge s derived states are found at the lower band edge. Lowering of the 
s states relative to the p states due to decreasing s-p hybridization across 
the series C-Si-Ge results in a depression of the optical band gap of GeN2 as 
compared to that of SiN2. Close coupling within the N2 pairs especially in 
CN2 gives rise to a third type of bonding, namely a- and 7r-type overlap of 
the N 2p states. The resulting bonding and antibonding bands are found at 
the upper edges of the valence and conduction bands of CN2 and again lead 
to a strong decrease of the optical band gap as compared to SiN2. 
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